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ABSTRACT 

A few Type la supernovae (SNe la) have been suggested to be an explosion of a 
super-Chandrasekhar-mass white dwarf (WD) to account for their large luminosi- 
ties, requiring a large amount of ^^Ni. However, the candidate over-luminous SNe la 
2003fg, 2006gz, and (moderately over-luminous) SN 1991T, have very different ob- 
servational features: the characteristic time-scale and velocity are very different. We 
examine if and how the diversity can be explained, by ID spherical radiation trans- 
port calculations covering a wide range of model parameters (e.g., WD mass). The 
observations of SN 2006gz are naturally explained by the super-Chandrasekhar-mass 
model. SN 1991T represents a marginal case, which may either be a Chandrasekhar or 
a super-Chandrasekhar-mass WD explosion. On the other hand, the low velocity and 
short time-scale seen in SN 2003fg indicate that the ejecta mass is smaller than the 
Chandrasekhar-mass, which is an apparent contradiction to the large luminosity. We 
suggest that the problem is solved if the progenitor WD, and thus the SN explosion, 
is aspherical. This may reflect a rapid rotation of the progenitor star, likely a con- 
sequence of the super-Chandrasekhar-mass WD progenitor. The observed differences 
between SNe 2003fg and 2006gz may be attributed to different viewing orientations. 

Key words: white dwarfs - radiative transfer - supernovae: individual (SN 2006gz, 
SN 2003fg, SN 1991T) 



1 INTRODUCTION 

Type la Supernovae (SNe la) are currently the most mature 
cosmological distance indicator which led to the discovery 
of the acceleration of the Universe (Riess et a. 1998; Perl- 
mutter et al. 1999). Their use as distance indicators relies on 
the well calibrated light curve characteristics, namely a phe- 
nomenological relation between the peak luminosity and the 
light curve width ("Phillips relation" or "stretching factor"; 
Phillips 1993; Perlmutter et al. 1997; Phillips et al. 1999). 

By clarifying natures of progenitor system(s) of SNe la 
(Livio 2000: Hillebrandt & Niemeyer 2000; Nomoto et al. 
2003), one expects to obtain deep knowledge on the origin 
of the light curve relation for the better luminosity calibra- 
tion, as well as new application of SNe la to a range of 
cosmological study. A progenitor of normal SNe la (Branch 
et al. 1993; Li et al. 2001) is believed to be a white dwarf 
(WD) having nearly the Chandrasekhar-mass. Recently, a 
few SNe la have been suggested to have originated from a 
WD whose mass exceeds the Chandrasekhar-mass of a non- 
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rotating WD, raising a possibility that not all SNe la are 
from a single class of progenitor system. 

The suggestion mainly relies on their exceptionally large 
luminosity, which requires more than IMq of ^^Ni synthe- 
sized and ejected during the explosion. The thermonuclear 
explosion converts the initial WD compositions partly to 
^''Ni. 7-rays from radioactive decay ^"^Ni ^^Co ^^Fe, 
through Compton scatterings and thermalization, power the 
optical emission of SNe la. Thus, the peak luminosity is 
closely related to the mass of ^®Ni initially synthesized (here- 
after MseNi). 

Historically, the first observationally-based suggestion 
for the super-Chandrasekhar WD explosion was made for 
SN la 1991T, as is summarized by Lira et al. (1998) and 
Fisher et al. (1999). However, the distance to this SN has 
been actively debated, and the smaller value is now favored 
(e.g., Saha et al. 2001). As a result, Msgni required to explain 
the peak bolometric luminosity has been reduced down; For 
the reddening E{B — y)host = 0.14 and the distance mod- 
ulus n = 30.74 as adopted by Stritzinger et al. (2006), the 
peak l^-band magnitude is Mv ~ —19.66. Stritzinger et al. 
(2006) derived the mass of '^''Ni, Akem ^ O.S7Mq, based on 
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the peak bolometric luminosity Lboi ~ 1-7 x 10 erg s^^. 
This amount of ^'^Ni could be marginally explained even 
by a Chandrasekhar-mass model, although one may need 
an extreme, purely detonation-driven explosion (Khokhlov, 
MMer, & Hoflich 1993). 

The progress had to wait until a discovery of over- 
luminous SNe la, which require that MseNi ^ IMq. The 
first candidate, SN 2003fg (SNLS-03D3bb), was reported by 
Howell et al. (2006). The SN reached Mv = -19.94 mag at 
peak and Lboi,pcak ~ 2.6 x lO'*'^ erg s"'^. From this, How- 
ell et al. (2006) estimated that MseNi ~ 1.29Mq, using 
an approximate relation between the peak luminosity and 
AfssNi- Since there should be other elements required both 
theoretically (i.e., all the well studied Chandrasekhar-mass 
models do not produce more than IMq of ^^Ni; Khokhlov 
et al. 1993; Iwamoto et al. 1999; Ropke et al. 2007 and 
references therein) and observationally (i.e., the maximum- 
light spectrum of SN 2003fg was dominated by intermedi- 
ate mass elements), Howell et al. (2006) argued that SN 
2003fg should come from a super-Chandrasekhar-mass WD. 
Two other over-luminous SNe la, thus candidates of the 
super-Chandrasekhar-mass WD explosions, have been re- 
ported: SNe la 2006gz (Hicken et al. 2007) and 2007if (Yuan 
et al. 2007). For SN 2006gz, Hicken et al. (2007) adopted 
E{B - V')host = 0.18 and fi = 34.95 (but also see §2), de- 
riving Mv — —19.74 at peak and Lbohpcak ^ 2.2 x 10*^ erg 
s-\ and thus M{^^m) ~ 1.2M0. 

In addition to the large luminosity, they are unique also 
in their spectra. Temporal sequence of the optical spectra of 
SN 2006gz is presented by Hicken et al. (2007). They iden- 
tified C II A4745, 5490, 6580, and 7324, in the spectra taken 
10—14 days before the maximum brightness, with the equiv- 
alent width reaching ~ 25A and at the velocity of ~ 15, 500 
km s~^. This is the strongest evidence for unburned carbon 
of any SNe la observed so far (SN 1990N by Leibundgut 
et al. 1991; SN 2006D by Thomas et al. 2007). Adding to 
this, a Si II 6355A absorption velocity was unusually low 
well before the maximum (~ 12, 500 km s^^ at ~ 10 days 
before the maximum). Without rapid decline observed in 
typical SNe la (Benetti et al. 2005), the Si II velocity at 
the maximum brightness settled down to usiii ~ 11, 500 km 
s~^, a value typical of normal SNe la. Based on this unusu- 
ally slow decline in the Si II velocity, Hicken et al. (2007) 
speculated that this might be a result of a dense unburned 
C-l-O envelope overlying a Si-rich region; The dense C-l-0 
envelope may be expected in the merger of two WDs lead- 
ing to the formation of the super-Chandrasekhar-mass WD 
progenitor. Howell et al. (2006) also reported the probable 
detection of a C II feature in the spectrum of SN 2003fg 
around the maximum brightness; if pre-maximum spectra 
had been available, the C II feature may have been in the 
spectra early on. These spectroscopic features suggest that 
the strong C II and the slow evolution of Si II velocity may 
be a distinct feature of over-luminous SNe la. However, the 
spectral sequence is only available for SN 2006gz, and it has 
not been clarified to what extent these features are common 
in over-luminous SNe la. Furthermore, it has not been clar- 
ified yet how these features, e.g., the dense C-l-O envelope, 
are related to the progenitor system (e.g., see Maeda et al. 
2008b for a caution of interpreting the C-|-0-rich region as 
an outcome of merging two WDs). 

Summarizing, the main argument for the super- 



Chandrasekhar model is the large peak luminosity and 
the mass of ^^Ni, with spectroscopic features only indica- 
tive. Detailed study on observed characteristics of super- 
Chandrasekhar-mass models is missing. Clearly, the lumi- 
nosity is not the only quantity that is dependent on the 
underlying models. The suggestion by Howell et al. (2006), 
about the Super-Chandrasekhar-model for SN 2003fg, is ac- 
tually based on another observed character, i.e., the velocity 
of the expanding SN materials as deduced from its spectra. 
SN 2003fg showed exceptionally small velocity around the 
maximum brightness (~ 8,000 km s~^) compare to that of 
normal SNe la (> 10,000 km s~^). At the first look, this 
seems fully consistent with the expectation, since the super- 
Chandrasekhar-mass WD provides the large binding energy, 
and thus the small kinetic energy per mass as compared to 
an explosion of a Chandrasekhar-mass WD. Jeffery, Branch, 
and Baron (2006) confirmed this statement on the relation 
between the WD mass and the velocity-scale. 

Although the argument may sound satisfactory, it is not 
a whole story. Studying emission processes should add fur- 
ther constraints on the underlying models. There are two 
over-luminous SNe la whose observed characteristics are 
available in the literature, i.e., 2003fg and 2006gz. Exten- 
sive observational data set for the moderately over-luminous 
SN 1991T is also available. By examining these observations 
(§2), we clarify that they have quite different properties. For 
example, SN 2006gz showed the velocity similar to normal 
SNe la - then, the argument for the super-Chandrasekhar- 
mass model for SN 2006gz may be flawedQ The aim of this 
paper is to examine if the observed characteristics can be ex- 
plained by super-Chandrasekhar-mass models, by means of 
radiation transfer calculations. We especially focus on SNe 
2003fg and 2006gz. In §3, we describe SN la models and 
summarize a method for the radiation transfer calculations. 
Results based on an extensive set of model calculations are 
shown in §4. Discussion is given in §5, where we examine un- 
certainties involved in our calculations. We also discuss the 
nature of these peculiar over-luminous SNe la. The emphasis 
is placed on our finding that not every observed character- 
istics can be interpreted within the context of the simple 
(spherical) super-Chandrasekhar-mass WD models, and on 
our suggestion that the progenitor should largely deviate 
from spherical symmetry to account for the observed fea- 
tures (at least for SN 2003fg). Implications on SN 1991T 
are also discussed. The paper is closed in §6 with conclud- 
ing remarks. 



2 OBSERVATIONAL CHARACTERISICS 

In this section, we summarize observed features relevant to 
this work, for SNe 1991T, 2003fg, and 2006gz. Howell et al. 
(2006) presented their observational data of SN la 2003fg. 
The multi-band light curves are fitted well by fc-corrected 
template light curves with the stretching factor (Perlmutter 
et al. 1997, 1999) of s = 1.13. The corresponding Amfg (the 
magnitude change in the first 15 days past B maximum: 
Phillips 1993) is the following: Amfg = 0.84 ± 0.02, using 

^ Indeed, it is shown in this paper that the low velocity seen in 
SN 2003fg has an apparent problem in the super-Chandrasekhar- 
mass model. 
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equation (5) of Perlmutter et al. (1997j^. We then convert 
Amfc, to t+i/2, the time since the maximum luminosity to 
half the maximum luminosity as measured in a bolometric 
light curve (Contardo, Leibundgut, & Vacca 2000). This is 
possible thanks to a correlation between these two quanti- 
ties. Note that i+i/2 in this paper is a post-maximum quan- 
tity, not the rise time frequently used in the related field. 
For the conversion, we fit an observationally derived set of 
(Amfg, t^i/2) for nearby 9 SNe la presented in Contardo et 
al. (2000), assuming a functional form of t+i/2 = aAmfs + fc. 
We obtain a = -4.52 ±0.786 and b = 17.3 ±1.03. Using this 
relation, t+i/a = 13.5 ± 1.7 days for SN 2003fg. The deriva- 
tion of t+i/2 in SN la 2006gz is straightforward. The UBVRI 
bolometric curve is presented by Hicken et al. (20070, and 
it directly gives i+i/2 = 18 days. The light curve width is 
very different between the two SNe la. 

They are very different in characteristic velocity, too. 
SN 2003fg showed a Si II absorption velocity (iisiii) of 
8,000 ± 500 km s"^ around maximum brightness (Howell 
et al. 2006). The same value for SN 2006gz is in the range 
11, 000 - 12, 000 km s'^ (Hicken et al. 2007). 

Typically, normal SNe la have f+i/2 in the range of 
~ 9 — 14 days (e.g., Contardo, et al. 2000) and vsiii in ~ 
10, 000-14, 000 km s"^ (e.g., Hachinger, Mazzah, & Benetti 
2006). If we focus on SNe la with MseNi ~ O.&Mq, the 
Phillips relation indicates that Amfg ~ 1.1 (Stritzinger et 
al. 2006). This corresponds to t+i/2 = 12.1 ± 1.8 using the 
relation we derived above. The examples are SNe 2003du 
and 1990N (e.g., Stritzinger et al. 2006), both having Usui ~ 
10,000 - 11,000 km s^^ 

The time scale, t+i/2, in SN 2003fg is at about the up- 
per bound of normal SNe la, while usiii is much lower, t+i/2 
in SN 2006gz is well above the range seen in normal SNe la, 
while usiii is similar to normal cases. As for the bolometric 
magnitude at the maximum light (I/boi,pcak), SN 2003fg has 
been claimed to be over-luminous; Lboi.pcak ~ (2.5 — 2.8) x 
10*^ erg s"^ (HoweU et al. 2006). SN 2006gz has also been 
claimed to be over-luminous; Z/boi.poak ~ (1.8— 2.6) x 10^^ erg 
s~^ assuming E{B — V')host ~ 0.18 and ^ = 34.95 (Hicken 
et al. 2007). However, there is a caveat for I/boi,peak in SN 
2006gz; Maeda et al. (2008b) pointed out that the peak 
luminosity of SN 2006gz derived by Hicken et al. (2007) 
involves large uncertainty, because the host-galaxy extinc- 
tion has not been well constrained. If the host extinction 
would be totally negligible, then Mv ~ —19.2 at peak and 
iboi.pcak ~ 1.3 X 10*'^ erg s"^ , which is relatively bright but 
not over-luminoufl 

For SN 1991T, the following values are relevant, t+i/2 = 
14.0 given by Contardo et al. (2000), based on compilation 
of muhi-band light curve, vsai ~ 10,000 - 11,000 km s"^ 
(Hachinger et al. 2006 and references therein) . In these prop- 
erties, SN 1991T represents an intermediate case between 

2 Hicken et al. (2007) derived Amfg ~ 0.9 for SN 2003fg in their 
Figure 4. Although this is different from what we derived here, 
the expected t+i/2 is ^ 13 days, which is well within the errors 
of our estimate 

^ The data file is presented at 

://www. cfa.harvard.edu/supernova/SNarchive.html I . 
* This uncertainty highlights the importance to investigate the 
radiation process and to obtain constraints independent from the 
peak luminosity, as is aimed in this paper. 



SNe 2003fg and 2006gz, and indeed close to those of normal 
SNe la. The peak luminosity with the favored distance mod- 
ulus is I/boi,poak ~ 1.7 X 10*'^ erg s~^ (Stritzinger et al. 2006) 
which is at the upper bound of normal SNe la, although the 
original suggestion for the super-Chandrasekhar-mass model 
was based on the brighter estimate (I/boi,poak ~ 2.3 x 10*^ 
erg s"^ Fisher et al. 1999; Contardo et al. 2000). 

Summarizing, the two most probable candidate super- 
Chandrasekhar-mass WD explosions, i.e., SNe 2003fg and 
2006gz, have very different observational characteristics, al- 
though this fact has not been emphasized in the literature. 



3 METHOD AND MODELS 
3.1 SN la Models 

We construct SN la models starting with the density struc- 
ture of the W7 model (Nomoto, Thielemann, & Yokoi 1984) 
which reproduces basic observational features of normal SNe 
la (Branch et al. 1985). Assuming the homologous expan- 
sion which should be a good approximation for an explosion 
of a compact progenitor, the density distribution as a func- 
tion of velocity is uniquely determined from the normalized 
reference density distribution (i.e., W7), by specifying the 
following model parameters: 

(i) Mwd: the mass of the WD. 

(ii) pc- the WD central density at the burning ignition. 

(iii) /ece: the mass fraction of Fe-peak elements pro- 
duced by strong electron captures, i.e., Fe-peak elements 
excluding ^''Ni (e.g., '^^Ni, ^'^Fe, '^''Fe). 

(iv) /ssNi: the mass fraction of ^^Ni. 

(v) /ime: the mass fraction of partially burned interme- 
diate mass elements like Mg, Si, and S. 

The mass fraction of unburned C±0 materials (/co) is then 
simply given by /co = 1 - /ece - /seNi - /ime. 

The procedure of the model construction follows that 
presented in Jeffery et al. (2006) (see also Howell et al. 2006 
and Maeda et al. 2008b). For given Mwd and pc, we com- 
pute the binding energy of the WD {Eh) by the formulae 
given by Yoon & Langer (2005), who examined a sequence 
of structure of a rotating WD. The energy produced by the 
nuclear burning (i?nuc) is given by Mwd, /ece, /seNi, and 
/ime, by a simple relation 

-B„uc/(10"'erg) = (1.74/ECE+1.56/56Ni±1.24/iME)Mwd/M0 .(1) 
The kinetic energy is then simply 

-Ek = Enuc — E\y . (2) 

The reference density structure is then scaled in a self-similar 
marmer, i.e., 

piv) oc Mli^E^'^' , (3) 
and 

^ « MJ^'e'/' , (4) 

where i; is a velocity of a Lagrangian fluid element, and p{v) 
the density there. 

Now that we have specified the density structure, we 
are left to specify the distribution of elements. We examine 
three extreme cases so that our models should cover the real 
situation anyway. 
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Figure 1. Examples of the density structure arc sliown for 
SW7IM (witli /ece = 0.1; solid) and for Sup2.3IM (with /ece = 
0.1; dashed). Also shown is the composition structure of SW7IM, 
for the three different mixing prescription. 



(i) Stratified (model sequence "a"): characteristic burn- 
ing layers are totally separated and stratified. The electron 
capture region (where the mass fraction of stable Fe-peak 
elements is set to be unity) is at the centre, surrounded by 
the ®^Ni-rich region (the mass fraction of ^®Ni is set to be 
unity), then by the partially burned layer (mass fractions 
of Si and S are set to be 0.7 and 0.3, respectively), and 
by the unburned C+0 layer at the outermost region (mass 
fractions of C and O arc set to be 0.5 and 0.5, respectively). 

(ii) Mixing in the Fe-rich region (model sequence "b"): 
Basically the same with the model sequence " a" , except that 
the innermost electron capture region and the ^^Ni region 
are assumed to be fully mixed in the composition structure. 

(iii) Fully mixing (model sequence "c"): The composi- 
tions are fully and homogeneously mixed throughout the 
ejecta. 

The models "a" and "b" are most likely the case at least in 
normal SNe la, in a viewpoint of observed spectral evolution 
(Stehle et al. 2005; Mazzah et al. 2008) . There are observa- 
tional hints favoring the models " a" than " b" for some SNe 
la (Hoflich et al. 2004; Motohara et al. 2006). Examples of 
the density and composition structures are shown in Figure 
1. 

We have examined an extensive set of models as shown 
in Table 1. Mwd is varied from 1.39Mq to 2.6Mq separ 
rated by O.SMq each, pc is set to be 3 x 10^ g cm~^ in 
most models, as the model sequence with varying pc can be 
self-similar (only different in -Ek). The dependence of the ob- 
served behaviors on pc is examined only for the models with 
Mwd = 2Mq, and can be generalized to different masses in 
a straightforward way. /seNi is fixed for given M^d, corre- 
sponding to MssNi = 0.6Mq for the Chandrasekhar-mass 
models (Af„d = 1.39Af©) and Ahem = I.OMq for the su- 
per Chandrasekhar-mass models (Mq > 1.39Mq). /ece 
is varied from (or 0.1 so that Ek should be positive) to 
1.0 — /56Ni, corresponding to the complete burning of whole 



ejecta to the Fe-peak elements, separated by 0.1 each. Ei- 
ther /iME or fco is set to be zero to reduce the number 
of model calculations. These two represent the two extreme 
cases, and thus should cover the real situation. The total 
number of models we construct here is then 246, including 
models with different /ece and difi^erent mixing prescrip- 
tion. 



3.2 Radiation Transfer 

For each model which is dcscritizcd in 172 radial zones, we 
performed radiation transfer calculations based on a numer- 
ical code presented in Iwamoto et al. (2000). The impor- 
tant assumption in our calculations is spherical symmetry. 
The transport of 7-rays from radioactive ®^Ni/Co/Fe de- 
cay is solved under a gray atmosphere approximation with 
— 0.027 cm'^ g~^, which should be a sufficiently accu- 
rate approximation for our present purpose (Maeda 2006a). 
With the energy input by 7-rays, the optical photon trans- 
port is solved again under the gray approximation. Opacity 
is provided by Thomson scatterings and line scatterings, i.e., 



(5) 



The ionization is determined by solving the Saha equation 
(i.e., in LTE), then used to determine fCe- at each time step 
and in each radial grid. 

The prescription of Kune is more complicated and uncer- 
tain, as it is affected by many weak lines whose atomic data 
arc experimentally not well known (Woosley et al. 2007). 
In this paper, we follow a simplified, phenomenological pre- 
scription based on modeling observations of normal SNe la, 
and restrict our selves to focus on bolometric light curves. 
The main advantage of this treatment is that the bolometric 
light curve is insensitive to a color evolution which can be 
quite complicated. We follow Mazzali et al. (2001): 

«:ii„e = 0.5x[0.25(Xece + X^wi) + 0.025(Xime + Xco)] cm^g"' .(6) 

Here, Xi denotes a local (position-dependent) mass fraction. 
The coefficient (0.5) is added to the original proscription 
of Mazzali et al. (2001) so that the SW7 model sequence 
satisfies the observational constraint of normal SNe la with 
MseNi ~ O.6M0 (§4). 



4 RESULT 

4.1 Fading Rate and Photospheric Velocity 

Figure 2 shows examples of the synthetic bolometric light 
curves. Shown in the figure are the models with the mixing 
"a" (stratified) and /ece = 0.1. If the other parameters axe 
similar, models with larger Mwd have broader light curve 
shape around the maximum brightness, i.e., t+i/2 is larger 
for larger Mwd. Also, the photospheric velocity at the max- 
imum luminosity (wpcak) is smaller for larger Mwd because 
of smaller -Ex/Mwd (note that the velocity in Supl.7IM is 
larger than SW7IM because of the larger amount of Mssni 
and thus larger -Ex/Mwd). 

This example shows that the relation between i+i/2 
and Vpeak can provide a strong tool to check the validity of 
a model as compared with the observations, although this 
has not been examined for SNe 2003fg and 2006gz in the 



Asphericity of Over-Luminous Type la Supernovae 5 



Table 1. SN la Model. 



Name 




Pc 


/ece 


/56Ni 


/iME 




fco 










SW7IM 


1.39 


3.0e9 


(0 - 0.57) 


0.43 


(0.57 - 


0) 







0.6 


1.40 - 


1.80 


LW7IM 


1.39 


3.0e9 


(0 - 0.28) 


0.72 


(0.28 - 


0) 







1.0 


1.53 - 


1.73 


Supl.yiM 


1.70 


3.0e9 


(0 - 0.41) 


0.59 


(0.41 - 


0) 







1.0 


1.55 - 


1.90 


Sup2IM 


2.00 


3.0e9 


(0 - 0.50) 


0.50 


(0.50 - 


0) 







1.0 


1.56 - 


2.06 


Sup2IM10 


2.00 


l.OelO 


(0 - 0.50) 


0.50 


(0.50 - 


0) 







1.0 


1.24 - 


1.74 


Sup2.3IM 


2.30 


3.0e9 


(0 - 0.57) 


0.43 


(0.57 - 


0) 







1.0 


1.56 - 


2.21 


Sup2.6IM 


2.60 


3.0e9 


(0 - 0.62) 


0.38 


(0.62 - 


0) 







1.0 


1.56 - 


2.36 


SW7CO 


1.39 


3.0e9 


(0 - 0.57) 


0.43 







(0.57 - 


0) 


0.6 


0.42 - 


1.80 


LW7CO 


1.39 


3.0e9 


(0 - 0.28) 


0.72 







(0.28 - 


0) 


1.0 


1.05 - 


1.73 


Supl.7CO 


1.70 


3.0e9 


(0 - 0.41) 


0.59 







(0.41 - 


0) 


1.0 


0.69 - 


1.90 


Sup2CO 


2.00 


3.0e9 


(0 - 0.50) 


0.50 







(0.50 - 


0) 


1.0 


0.32 - 


2.06 


Sup2CO10 


2.00 


l.OelO 


(0 - 0.50) 


0.50 







(0.50 - 


0) 


1.0 


0.001 


- 1.74 


Sup2.3CO 


2.30 


3.0e9 


(0.10 - 0.57) 


0.43 







(0.47 - 


0) 


1.0 


0.33 - 


2.21 


Sup2.6CO 


2.60 


3.0e9 


(0.10 - 0.62) 


0.38 







(0.52 - 


0) 


1.0 


0.01 - 


2.36 



Models with the W7 reference density distribution. For the description of the parameters, see the main text. Each column represents a 
set of models: For example, model SW7IM is examined with three different mixing prescription (mixing "a", "b", and "c"; see the main 
text). For each mixing, a set of (/ece, Iime) - (0, 0.57), (0.1, 0.47), (0.2, 0.37), (0.3, 0.27), (0.4, 0.17), (0.5, 0.07), and (0.57, 0) - are 

examined. 
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Figure 2. Examples of the synthesized bolometric light curves. 
Shown here are the models with the mixing "a" (stratified) 
and /ece = 0.1; SW7IM (thick-solid), LW7IM (thin-solid), 
Supl.7IM (thick-dashed), Sup2.0IM (thin-dashed), Sup2.3IM 
(thick-dotted), and Sup2.6IM (thin-dotted). Shown in parenthe- 
ses are the photospheric velocity at the maximum bolometric 
brightness (fpcak km s~-^). 

previous studies. The constraints are the following: (1) A 
model should reproduce derived in §2. (2) The pho- 

tospheric velocity at the maximum luminosity in a model 
(Vpeak) should be equal to or less than observed wsiii, since 
the Si II absorption must be formed above the photosphere. 

Figure 3 shows (t+i/2, Wpoak) for all the models, ex- 
cept for Sup2.0IM10 and Sup2.0COlcE]. The SW7IM/CO 



5 We have found that (t+i/2, -Upoak) of Sup2.0IM10/Sup2.0CO10 
can be obtained by simply shifting those of Sup2.0IM/CO, along 



model sequence (with AfseNi = O.GA/q) approximately sat- 
isfies the (i+i/2, I'Siii) constraint of the normal SNe la with 
MseNi ~ O.6M0 (e.g., SNe 2003du, 199QN, as the observed 
range shown in Figure 3) . 

Important discovery in Figure 3 is that even a sin- 
gle model does not satisfy the observed characteristics of 
SN 2003fg. On the other hand, SN 2006gz can be ex- 
plained by many models (note that observed iisiii can be 
larger than synthetic Wpcak), and most naturally by Super- 
Chandrasekhar models having Msbni ~ IA/q, which results 
in iipeak ~ t'Siii for the favored non-mixing case ("a"). We 
have examined an extensive set of models including extreme 
cases, so that our failure to reproduce the observed charac- 
teristics of SN 2003fg represents a principal difficulty in the 
super-Chandrasekhar model for this SN la. 

We found that the lower boundary in the t^i/2 — t^pcak 
plot, covered by varying various parameters for given Afwd, 
is represented by a model sequence of " a" with only varying 
/ece- In the mixing cases ("b" and "c"), ^®Ni is mixed down 
to the low velocity, and thus the diffusion time-scale becomes 
long as compared with the stratified case "a". This results 
in larger t+i/2 in mixing "b", "c" than "a". We also found 
that the curves, obtained by varying /ece, almost overlap 
between the cases "IM" (/co = 0) and "CO" (/ime = 0) in 
the t+i/2 — Wpeak plane, and the ones with /ime = (e.g., 
SW7CO) cover the wider range of parameter space. 

Thus, the model sequence with the mixing "a" and 
fiME = form the lower boundary as we vary the value for 
/ece, in the t-t-i/2 — Wpcak plane for given Afwd and AfseNi. 
If the observationally derived set of (f+i/2, vsni) is above 
the curve, the model sequence is acceptable (i.e., there are 
observationally acceptable combinations of parameters for 
given Afwd and AfssNi). If it is below the curve, the combi- 
nation of Afwd and AfssNi should be rejected. Figure 4 shows 
the lower boundary of Upcak as a function of t+i/2 for dif- 

the curve defined by the mixing "a" to the bottom-right direction. 
This is because larger pc is equivalent to lower Ej<^ . Thus changing 
Pc can be mimicked by changing (i-e., by changing /ece)- 



6 K. Maeda and K. Iwamoto 



20000 
18000 
16000 
t7 14000 

CO 

E 

12000 
>" 10000 
8000 
6000 



Normal 




. f,':S AO C'^* A 



o oft 



SN 199lt ' ° 



SN 2003fg SN 2006gz 



10 



20 
t.i. [day] 



30 



20000 
18000 
16000 
14000 
12000 
10000 
8000 
6000 




Normal 



SN 2003fg 



10 



20 



30 



t [day] 



Figure 3. {t+i/2i fpcak) shown for the models. The meaning 
of the symbols are the following: SW7IM and SW7CO (filled 
squares), LW7IM and LW7CO (open squares), Supl.7IM and 
Supl.7CO (filled circles), Sup2.0IM and Sup2.0CO (open cir- 
cles), Sup2.3IM and Sup2.3CO (filled triangles), Sup2.6IM and 
Sup2.6CO (open triangles). The SW7 models with different /ece 
are connected by the lines, separately for different mixing pre- 
scription ("a", "b", and "c"); solid lines for SW7IM and dotted 
lines for SW7C0. The observed range for normal SNe la with 
A^56Ni ~ O.6M0 is shown by the shaded region. The observed 
characteristics of SNe 2003fg, 2006gz, and 1991T, are also plot- 
ted as the shaded regions with observed uncertainty evaluated in 
§2. 

ferent Mwd (note that Msbni is fixed for given Af„d in this 
paper) . 

Figure 4 shows that the curve is similar for models hav- 
ing the same Msbni- The behavior can be explained as fol- 
lows to the first approximation. The time-scale of the hght 
curve evolution around peak is scaled as follows (Arnett 
1982). 

, 1/2 « f3/4 ^-1/4 ,„\ 

where k is the opacity averaged over the ejecta. The estimate 
of the photospheric velocity at the maximum brightness is 
complicated. For the present demonstration purpose, it is 
enough to assume that this is scaled with the average veloc- 
ity, i.e., 

^pcak « M-I'^EU^ . (8) 

Then, combining these two expressions, we obtain the fol- 
lowing relation: 

fpeak (X + • (9) 

The opacity is provided by Fe-peak elements, and stable 
Fe-peak elements below the ^®Ni-rich region does not sig- 
nificantly contribute to the opacity in non- mixing case "a". 
Thus, 

llpcak tX M56Nit+i/2 ' (10) 

Therefore, the lower boundary (i.e., non-mixing case "a") 
in the — tipcak plane for a given model is basically 



Figure 4. The lower boundary in the t+i/2 — ''peak plane shown 
for the W7-sequence models with various M^^. The meaning of 
the symbols are the following: SW7CO (filled squares and thick 
solid line), LW7CO (open squares and thin solid line), Supl.7CO 
(filled circles and thick dashed line), Sup2.0CO (open circles and 
thin dashed line), Sup2.3CO (filled triangles and thick dotted 
line), Sup2.6CO (open triangles and thin dotted line). For the 
observed values (shaded regions), see the caption of Figure 3. 



determined by MseNi. Note that the curves representing 
the lower boundary in Figure 4 do not exactly behaves 
as Upcak oc t^y2' because the models with different /ece 
are not self-similar to one another. Specifically, a model 
with lower /ece (thus larger t+i/2 and lower «peak) has 
i-i-i/2 larger than expected by this simple analytic relation, 
since they have more centrally concentrated ^®Ni distribu- 
tion leading to the larger diffusion time-scale. 

One may wonder if our results are consistent with obser- 
vations of normal SNe la and some peculiar SNe la. First of 
all, we emphasize that we restrict ourselves to investigate 
Chandrasekhar models with Msbni O.GAf© and super- 
Chandrasekhar models with Msbni ^ lAf© . As such, our cal- 
culations do not directly cover (a) detailed variation in nor- 
mal SNe la, and especially (b) sub-luminous SNe la which 
requires Msbni < 0.6Mq. 

However, these issues are consistent with our results, 
and thus it is safe to use our calculations to model the super- 
Chandrasekhar SNe la candidates: (a) A range of t+i/2 seen 
in normal SNe la can be explained by varying MseNi. For 
example, if we take Msbni < 0.6Mq, the model sequence 
moves to the left [i.e., to smaller i+i/2; equation (10)]. This 
also reduces the peak luminosity as ^^Ni is the energy source. 
In this way, (t+i/2, fsiii) of normal SNe la can basically 
be explained, simultaneously satisfying the Phillips relation 
(Mazzali et al. 2001). (b) There are two SNe la showing 
Wpcak smaller than SN 2003fg. They are "sub-luminous" pe- 
culiar SNe la 2002CX (Li et al. 2003; Branch et al. 2004; 
Jha et al. 2006) and 2005hk (Jha et al. 2006; Phillips et 
al. 2007; Sahu et al. 2008), suggested to be a less ener- 
getic explosion with smaller mass of ^^Ni than normal. SN 
2005hk showed usiii ~ 6,000 km s~^ i^pcak according to 
the spectrum modeling; Sahu et al. 2008) and t+i/2 ^ 20 
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days. Some of our models have small Ek (Tab. 1), which 
results in small iipcak and large t+i/2- For example, Figure 4 
shows that for sufficiently small Ek, Wpeak ~ 8,000 km s~^ 
and t+i/2 ~ 30 days for A/sbni ~ IMq. It has been esti- 
mated that MseNi ~ 0.2Mq in these faint SNe la (Phillips 
et al. 2007; Sahu et al. 2008), and a model with such a 
small amount of ^^Ni should be considered. From equation 
(10), we see that t+i/2 c>c M^^^^v'^^^ . Applying our result 
to this relation, we expect that i+1/2 ~ 16 — 20 days for 
AfssNi ~ 0.2 — 0.3A/q and for Wpcak ~ 6,000 km s^^, as is 
consistent with the observational features of SN 2005hk. 

Equations (9) and (10) are only illustrative, but already 
provide a solid argument which was confirmed by a set of 
our model calculations: (1) For given Afwd, if Ek is larger, 
then t+i/2 is smaller and Wpeak is larger. (2) To reproduce 
the observationally derived value of f+i/2, one can thus play 
with Ek (in terms of /ece, for example). The model at the 
same time predicts smaller Upcak for smaller MsgNi- 

The combination (t+i/2, vsm) of SN 2006gz is totally 
consistent with the expectation from MssNi ~ IMq. On 
the other hand, the above relation between t+i/2 and i^poak 
raises a difficulty in interpreting the observational data of SN 
2003fg, since (t+i/2, wsiii) falls into the range even below the 
lower boundary of the Chandrasekhar model with MsgNi = 
0.6Mq. 

4.2 Peak Luminosity and Phillips Relation 

The t+i/2 — iipeak relation is the main point in the present 
paper. We also examine the bolometric peak luminosity 
(^boi,peak) for the self-consisteucy test in this section. 

Figure 5 shows Lboi.poak as a function of t+i/2 for all 
the models (except for Sup2.0IM10/CO10). Normal SNe la 
(shown in stars in Fig. 5) can basically explained by the SW7 
model sequence, if we allow dispersion in Msbni (i.e., I/poak oc 
AfsBNi if other properties are the same). Our model value 
(MseNi ~ IMq) is slightly smaller than required by the peak 
luminosity of SNe 2003fg. The best value is MseNi ~ 1.2Mq, 
as is consistent with the estimate by Howell et al. (2006). For 
SN 2006gz, the best value is MseNi ~ I.IMq, if we adopt 
E{B - V)host = 0.18 and Rv = 3.1 (Hicken et al. 2007). 
The stringent lower limit for SN 2006gz is given by setting 
the host extinction negligible (§2); MseNi ^ 0.7Mq. 

With these values, let us go back to the discussion on 
the t+i/2 — «peak relation. If we take Msbni ~ I.IMq for SN 
2006gz, the lower boundary in Fig. 4 moves slightly to the 
right, but still marginally consistent with SN 2006gz. For SN 
2003fg, further increasing A/ssni makes the situation worse: 
the lower limit of Upcak now increases, making the deviation 
larger. Thus, the conclusion in the previous section does not 
change. 

So far we have seen that properties of SN 2006gz is 
consistent with super-Chandrasekhar models. Then, can we 
identify a set of model parameters (Mwd and MseNi) rele- 
vant to this SN? The t+i/2 — ^Siii constraint can be sat- 
isfied as long as Msbni ^ I.IMq. This is only the upper 
limit corresponding to the condition Upcak ^Siii- On the 
other hand, Lpoak requires MseNi -SJ IMq for Av ~ 0.56 
and MssNi ^ 0.7Mq for Av = 0. Thus, we conclude that 
an explosion producing Msbni ~ 1 — I.IMq is the most 
likely explanation for SN 2006gz, with possible range of 
MseNi ~ 0.7 — IAAIq. Directly deriving Afwd, rather than 
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Figure 5. ibol.poak ^ function of t^x/2- See the caption of 
Fig. 3 for the meaning of the symbols. For SN 2006gz, shown 
here are the range of the luminosity from Hicken et al. (2007) 
with likely Ay = 0.56 for the host extinction, and the strin- 
gent lower limit with Av = (see Maeda et al. 2008b). SN la 
samples from Contardo et al. (2000) are shown by open stars 
(among which SN la 1991T is explicitly indicated, but with the 
luminosity from Stritzinger et al. 2006). The observationally- 
derived relation between the t-|-i/2 and Lbol.poak (equivalent to 
the Phillips relation) is shown by a solid line (excluding bright SN 
la 1991T) and by a dashed line (excluding bright SN la 1991T and 
faint SN la 1991bg). An approximate linear fit to the results of 
the Sup2.0 sequence (with mixing "a") is indicated by the solid 
line. Two lines ("xl.l" shown by dashed line, "xO.8" shown by 
dotted line) roughly indicate what is expected for models with 
A^56Ni = 1-1A^0 and O.SMq, respectively. 

MseNi, turns out to be difficult: Af„d is degenerated in the 
i+i/2 ^ fpcak relation. In terms of Lpoak, Mwd is again ba- 
sically degenerated according to the approximate relation 
£'boi,pcak oc MseNi • Although there is a diversity resulting 
from different Mwd, the model luminosity differs up to only 
~ 10 per cent for the same mixing prescription, between 
Mwd ~ 1.4Mq and 2.6Mq. The variation is smaller than 
the observed error in Lboi,pcak (~ 20 per cent), and thus it 
is difficult to derive Mwd from the analysis presented in this 
paper. 

4.3 Temperature and Peak Luminosity 

Additional test on the models can be provided by the ejecta 
temperature. The discussion in this section is only qualita- 
tive: deriving the photospheric temperature requires detailed 
spectrum modeling (e.g., Hachinger et al. 2008), which is be- 
yond the scope of the present study. 

The ratio of equivalent widths of Si II A5972 to Si II 
A6355 provides a temperature indicator. In normal SNe la, 
the ratio is smaller for SNe with larger peak luminosity, 
and thus higher temperature (Nugent et al. 1995; Hachinger 
et al. 2008). The ratio in SN 2003fg is similar to that in 
relatively faint SN la 1994D (fig. 3 of Howell et al. 2006), 
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Figure 6. The effective temperature at bolometric maximum 
briglitncss, as a function of The meaning of the symbols 

are same as in Figure 4. Models shown here are for the mixing 
case "a". 



and that in SN 2006gz is similar to normal SN la 2003du (fig. 
1 of Hicken et al. 2007). These infer that the temperature of 
SN 2003fg is lower than normal SNe la with Msbni ~ O.GMq, 
and that in SN 2006gz is comparable to the normal case. 

Figure 6 shows the effective temperature at bolometric 
maximum brightness in our models. The effective temper- 
ature is generally a decreasing function of i+i/2; Although 
the photospheric velocity is smaller for larger i+i/2 (Fig. 4), 
the peak luminosity is smaller and the photospheric radius 
(i.e., the photospheric velocity multiplied by the peak date) 
is larger for larger t+i/2- The effect of the latter functions is 
more important, resulting in the dependence in Figure 6. For 
models resulting in same t+i/2, the temperature is larger for 
larger Mwd (see Fig. 4: The photospheric velocity is slightly 
smaller for larger M„d)- 

The temperature of SN 2003fg looks to be lower than 
normal ( i.e., SW7 models), and i+i/2 is comparable to nor- 



mal SNe la (i.e., t 



+1/2 ~ 12 — 14 days). This requires that 
the ejecta mass is smaller than in SW7. This is again con- 
sistent with our conclusion independently derived using the 
photospheric velocity. The temperature of SN 2006gz is com- 
parable to normal SNe la, but t+i/2 (~ 18 days) is larger 
than normal cases. This indicates that the ejecta mass is 
larger than in normal SNe la (Fig. 6). 



5 DISCUSSION 

5.1 Evaluation of Uncertainties 

One may wonder that uncertainty involved in our calcula- 
tions might change our conclusion. In this section, we exam- 
ine dominant sources of the uncertainty in our calculations. 
Here we examine effects of (a) different opacity prescription 
and (b) density distribution. In short, these do not alter our 
conclusion. 



Figure 7. The same as Figure 4, but for different opacity pre- 
scription, independent from composition. 



5.1.1 Opacity 

We have performed the same calculations for all the mod- 
els, with different opacity prescription. Here, the opacity 
takes the following form: Km-ic = 0.05 cm'^ g^^, indepen- 
dent from composition. The result is shown in Figure 7. 
The combination {t^i/2, tipeak) of the Super-Chandrasekhar 
WD model is again above that of the Chandrasekhar model 
and above (t+i/2, wsiii) of SN 2003fg. The observed val- 
ues of SN 2006gz are slightly above the lower boundary of 
the Super-Chandrasekhar model, and thus this SN is quite 
consistent with the Super-Chandrasekhar explosion. In this 
case, Upcak oc Mwd as the opacity is independent from the 
composition. The same conclusion as in the previous section 
applies, only if we use Mwd as the major function rather 
than MsBNi: To explain the observational characteristics of 
SN 2003fg, we require Mwd smaller than in the normal SNe 
la. 



5.1.2 Density Distribution 

The models with exponential or constant density distribu- 
tions are examined for selected model parameters to check 
the uncertainty. Most of existing explosion models (includ- 
ing the W7 model) predict more or less exponential density 
distribution (e.g., Woosley et al. 2007). We constructed the 
exponential density distribution from Mwd and Ek accord- 
ing to description in Jeffery et al. (2006) and Woosley et al. 
(2007). The construction of the constant density distribu- 
tion is trivial. These two types of density distribution are 
examined only for the following cases: 

(i) Mwd = 1.39M0 and 2M0 

(ii) pc = 3 X 10^ g cm"^ 

(iu) /56Ni 

2M^ 



0.43 (Mwd = 1.39Mq) or 0.50 (Mwd 



^0. 
(iv) 



/iME = 

mixing 
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Figure 8. The same as Figure 4, but for models with expo- 
nential or constant density distribution. Shown here are the 
following: exponential distribution with M„(j = 1.39Mq and 
-'^56Ni = 0.6Mq (filled circles and thick dashed line), exponen- 
tial with M^(j = 2Mq and MsgNi = (open circles and 
thin dashed line), constant density with Af^d = I-SQMq and 
-'^56Ni = O.6M0 (filled triangles and thick dotted line), con- 
stant with M„(j = 2Mq and MsgNi = (open triangles and 
thin dotted line). For comparison, the W7-sequence is also shown 
(SW7CO by filled squares and thick solid, Sup2.0CO by open 
squares and thin solid). 

For each model, /ece is varied to form the lower boundary 
in the f+1/2 - '^pcak plot for given M„d. 

The result is shown in Figure 8. Exponential models 
are similar to the W7 density models in the — Vpak 

plane. The constant density distribution predicts smaller 
fpeak than in the other models, but not drastically change 
the result. The constant density distribution is an extreme 
assumption, probably resulting in the lowest value of i^poak 
for given Thus, the uncertainty in the density distri- 

bution, as long as spherical symmetry is assumed, does not 
change our conclusion. 

5.2 SN 2003fg 

In this section (§5.2) and the following two sections (§§5.3 
& 5.4), we summarize our results for individual objects, and 
discuss implications. For SN 2003fg, we have found that the 
observational characteristics cannot be put into the super- 
Chandrasekhar-mass WD explosion scenario; Most impor- 
tantly, the observed — t^pcak relation (§4.1) requires that 
either Msbni or M^^ (or both) should be smaller than even 
the Chandrasekhar mass, contrary to the earlier expecta- 
tions (Howell et al. 2006). Additional support is provided by 
the ejecta temperature (§4.3), which also indicates that the 
ejecta mass (i.e., Mwd) is smaller than the Chandrasekhar 
mass. On the other hand, the large peak luminosity requires 
that Mwd ~ I.IMq (§4.2; Howell et al. 2006). 

This is apparently a contradiction. To remedy the prob- 
lem, we suggest that the ejecta structure is far from spheri- 
cal. Our radiation transfer calculations in this paper assume 
spherical symmetry, and for example proportional coefficient 



in equation (8) for «pcak should be a function of the view- 
ing angle in the presence of large deviation from spherical 
symmetry. If the viewing angle is such that the effective 
(isotropic) mass is small along the line-of-sight, this may ef- 
fectively look like an explosion with small amount of ^^Ni 
and/or Mwd, although the large luminosity can probably be 
provided by ^®Ni in the whole ejecta, not only along the 
light-of-sight (note that the photospheric velocity is more 
sensitively affected than the peak luminosity in asymmetric 
SN models; Maeda et al. 2006b; Tanaka et al. 2007). 

The important finding here is that the mass should be 
effectively small as viewed from an observer to satisfy the 
t_i_iy2 — fpeak Constraint. Such an explosion can not be a 
strongly jetted explosion of a spherically symmetric progen- 
itor star: The jet-type explosion should yield smaller t+i/2 
for an observer closer to the jet direction (thanks to the large 
isotropic Ek/Mwb) but at the same time resulting in larger 
^poak (for the same reason) (see e.g., Maeda et al. 2006b; but 
see Hillebrandt, Sim, & Ropke 2007). Alternatively, we sug- 
gest that the progenitor star is highly aspherical. This may 
actually be consistent with the Super-Chandrasekhar model, 
as such a massive WD should rotate rapidly to support an 
exceeding mass. 

Although an explosion based on such a deformed con- 
figuration has not been examined except for the purely det- 
onation model (Steinmetz, Miiller, & Hillebrandt 1992), we 
believe it is rational to assume that the structure after the 
explosion preserves the initial configuration to some extent. 
The disc-like structure has effectively small isotropic mass if 
viewed along the axis of the rotational symmetry. We sug- 
gest this is a situation in SN 2003fg. 

Hillebrandt et al. (2007) suggested an off-centre explo- 
sion model (a kind of one-sided jet-like explosion model) 
within the context of a Chandrasekhar-mass WD explosion, 
as an alternative explanation for SN 2003fg (see also Sim 
et al. 2007). Their argument is based on the viewing angle 
effect on the light curve features, i.e., the larger luminos- 
ity and smaller diffusion time-scale for an observer closer to 
the direction of the ^^Ni-rich blob. We emphasize the dif- 
ferences between our and their works: (1) Our suggestion 
on the ejecta asymmetry is based on the combination of the 
light curve and spectral features, and (2) the disc-like/oblate 
geometry, as we favor in this work, is different from their 
suggestion. 

Although we have clarified the need for the ejecta asym- 
metry for SN 2003fg, the progenitor WD mass (Mwd) is 
not conclusively constrained by the present study. This will 
need the detailed multi-dimensional radiation transfer cal- 
culations; especially, we need to understand the dependence 
of the luminosity on the geometry and the viewing angle. 

5.3 SN 2006gz 

We have found that observed features of SN 2006gz are 
consistent with expectations from the super-Chandrasekhar- 
mass WD explosion scenario. A difficulty in identifying the 
underlying model for SN 2006gz is the uncertainty of its 
host's extinction and the luminosity. Although we cannot 
perfectly reject the possibility that this is an explosion of a 
Chandrasekhar WD, we favor that this is indeed a super- 
Chandrasekhar WD explosion; the f+i/2 — Wpoak relation 
(§4.1) and the ejecta temperature (§4.3) are consistent with 
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the expectations from a super-Chandrasekhar WD explo- 
sion. 

For SN 2006gz, it is not necessary to introduce the 
ejecta asymmetry, unlike SN 2003fg. These two SNe may be 
intrinsically different, or, SN 2006gz may be explained by the 
configuration similar to that of SN 2003fg, but viewed at a 
different angle. Here, we point out that the difference in the 
observational features of these two SNe are qualitatively con- 
sistent with expected difference arising from the same con- 
figuration but viewed at different orientations. First, these 
two SNe la require similar amount of ^''Ni as estimated from 
the peak luminosity, and the main difference is in the spec- 
troscopic feature (upoak) and the light curve width (t+i/2)- 
The latter two are expected to be more sensitively dependent 
on the viewing angle than the estimated value of M(^®Ni) 
(Maeda et al. 2006b). For the disc-like/oblate configuration, 
we expect larger Wpcak and larger for the larger incli- 

nation. These two values in SNe 2003fg and 2006gz follow 
this tendency. This is, however, only a qualitative argument, 
and discriminating these two possibilities (the difference in 
the ejecta shape or the viewing angle) needs more detailed, 
multi-dimensional study. 

Maeda et al. (2008b) reported that SN 2006gz is pe- 
culiar also at late phases, > 300 days after the explosion. 
Applying the standard ^^Ni/Co/Fe heating scenario, they 
estimated MseNi by a factor of 5 smaller that estimated 
with the early-phase dat4f|. This contradiction could mean 
either of the following two possibilities; (f) SN 2006gz was 
not powered by decays of ^^Ni/Co/Fe at the early-phases 
(which then casts doubt on the super-Chandrasekhar-mass 
WD progenitor), or (2) it was powered by the decay, but for 
some reason (e.g., thermal catastrophe or dust formation), 
the bulk of the emission might be shifted to NIR/Mid-IR. 
Based on the result in this paper, we favor the second pos- 
sibihty. 



5.4 SN 1991T 

For SN i99iT, t+i/2 is at the upper boundary of normal 
SNe la, and usiii is at the lower boundary. This is marginally 
explained by the SW7 model sequence with Msbni = O.GMq . 
Thus, this constraint results in A/sbni ^ 0.6Mq. On the 
other hand, the other constraint from I/boi,pcak results in 
MseNi > O.SMq (Fig. 5). 

Strictly speaking, these two constraints are not mutu- 
ally consistent, and the same argument for SN 2003fg can 
apply to SN i99iT. However, the deviation between the ob- 
servations and models is not as large as in SN 2003fg, and 
SN f 99f T may be marginally consistent with a spherical ex- 
plosion with MseNi ~ O.SMq within uncertainties involved 
in our model calculations. 

Alternatively, SN f 99iT may also be an aspherical, disc- 
like explosion like our suggestion for SN 2003fg. However, the 
intrinsic property (e.g., the explosion geometry, the mass 
of ^"^Ni) of SN i99iT is likely different from that of SN 
2003fg (note that there is a possibility that SN 2006gz is 
intrinsically similar to SN 2003fg). The estimated value of 
M(5^Ni) is different, and other features (wpcak, t+1/2) do 



This is independent from the uncertainty in the host extinction. 



not seem to be consistent with the expectation from the 
disc-like/oblate geometry which we favor for SN 2003fg. 



6 CONCLUDING REMARKS 

In this paper, we critically examined super-Chandrasekhar- 
mass WD models for (candidate) over-luminous SNe la 
2003fg, 2006gz, and moderately over-luminous SN fa i99f T. 
Our new approach is to use two observed features, and 
'^pcak. This is equivalent to make use of a combined informa- 
tion of the light curve and spectra, and thus is more powerful 
than previous studies which mainly relied on the peak lumi- 
nosity. Our conclusions are summarized as follows: 

(i) Somewhat negatively, the observations of SN 2003fg 
are not readily explained by the standard ^^Ni/Co/Fe heat- 
ing scenario. The combination of relatively small t+i/2 and 
small «pcak requires that either Mwd or MseNi (or both) 
should be smaller than in normal SNe la. 

(ii) The observations of SN 2006gz can be naturally ac- 
counted for, by the (spherical) super-Chandrasekhar-mass 
model with MseNi ~ f — f.fAfQ. Although the peak lumi- 
nosity of SN 2006gz is largely uncertain (§2), we found that 
the super-Chandrasekhar-mass model is consistent with var- 
ious observational features of SN 2006gz, and thus we favor 
the interpretation that SN 2006gz was indeed over-luminous 
at peak as suggested by Hicken et al. (2007). 

(iii) Observed features of SN f99fT is marginally ex- 
plained by a spherical explosion of a WD and M(^^Ni) 
~ O.8M0. This may either be a Chandrasekhar WD or a 
super-Chandrasekhar explosion. 

The failure of fitting SN 2003fg is not a result of the 
M„d — Ek relation expected for SNe fa. In short, the large 
amount of ^''Ni should inevitably yield too large diffusion 
time-scale. One can then try a model with large Ek to re- 
duce the diffusion time-scale, but this inevitably leads to 
the large velocity - thus the contradiction. Therefore, the 
observed characteristics of SN 2003fg are indeed inconsis- 
tent with any parameter set within the standard ^''Ni/Co/Fe 
heating scenario. 

We suggest a solution to remedy the problem - the 
ejecta asymmetry resulting from a disc-like or oblate density 
structure of a rapidly rotating progenitor WD. This may in- 
deed be consistent with the super-Chandrasekhar WD sce- 
nario. Within the results obtained by the present calcula- 
tions, however, we can not conclusively derive the progenitor 
WD mass. An interesting possibility is that the different ob- 
servational features of SNe 2003fg and 2006gz may be unified 
into one scheme, i.e., an explosion of super-Chandrasekhar 
WD viewed at different directions. 

Thomas et al. (2002) pointed out that a signature 
of ejecta asymmetry can be identified in the absorption 
strength of Si H A6355. They concluded that the typical 
scale of inhomogeneity should be smaller than the size of 
the photosphere in normal SNe la, in order to account for 
uniformity of the absorption depth of the Si H. From fig. f of 
Hicken et al. (2007), we see that the absorption strength is 
indeed different between SNe 2003fg and 2006gz; SN 2003fg 
seems to show the weaker absorption. This may support the 
hypothetical configuration mentioned above, i.e., disc- like 
ejecta structure with SN 2003fg viewed closer to the polar 
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than SN 2006gz; In the polar direction, there are a small 
amount of materials to absorb the light emitted at the pho- 
tosphere, leading possibly to the weaker absorption strength 
in SN 2003fg. 

Our results can be tested by future observations of over- 
luminous SNe la. First, we expect a large degree of polariza- 
tion in SN 2003fg-like over-luminous SNe la, depending on 
the viewing direction. If SN 2006gz is an explosion similar 
to SN 2003fg, but viewed at different orientation, we also 
expect a large degree of polarization for SN 2006gz-like SNe 
la. Next, late-time nebular spectra may be used to infer the 
distribution of ^'^Ni directljQ (Motohara et al. 2006; Maeda 
et al. 2008a; Modjaz et al. 2008). Furthermore, if we plot a 
number of over-luminous SNe la on the f+i/2 — '^pcak plot, 
it will tell in a statistical way how large the deviation from 
spherical symmetry on average is, and how much the intrin- 
sic diversity (e.g., A/„d) in the super- Chandrasekhar-mass 
WD explosions is. 

More detailed, multi-dimensional studies are necessary 
to confirm our suggestions and speculations. First, multi- 
dimensional hydrodynamic simulations, based on a rapidly 
rotating super-Chandrasekhar WD, should be useful to 
check if our preferred geometry, i.e., the disc-like or oblate 
distribution of density and ^''Ni, can indeed be realized. 
Next, multi-D radiation transport calculations should be 
useful to see if the t+i/2 — Wpcak constraint can indeed be 
satisfied. We postpone these studies to the future, since in 
this paper we concentrate on examining a large parameter 
space to clarify the applicability and difficulty in spherically 
symmetric models (practically impossible to do in expen- 
sive multi-dimensional study), and to clarify the need for 
the non-spherical ejecta (for SN 2003fg). 

Identifying the geometry should be useful to evaluate 
possible contamination of SNe la, which do not follow the 
Phillips relation, in the cosmological study. SN 2006gz is 
roughly consistent with the Phillips relation if we adopt the 
peak luminosity derived by Hicken et al. (2007) (§4.2 and 
§4.3; Fig. 5), while SN 2003fg is not (Howell et al. 2006). 
Our scenario for the disc-like/oblate ejecta infers that over- 
luminous SNe la which do not follow the Phillips relation 
(e.g., SN 2003fg), should also show the peculiar low photo- 
spheric velocity. Thus, by performing spectroscopy, the SNe 
la with large deviation from the Phillips relation can be 
easily identified, and can be removed from samples for the 
cosmological study. 
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